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Measurements of the properties of Ac(2595), Ac(2625), Sc(2455), and Sc(2520) baryons 
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We report measurements of the resonance properties of Ac(2595)''' and Ac(2625)'*' baryons in their 
decays to A^7r"^7r~ as well as Ec(2455)^''''° and Ec(2520)'^^''' baryons in their decays to A^tt* 
final states. These measurements are performed using data corresponding to 5.2 fb~^ of integrated 
luminosity from pp collisions at ^/s — 1.96 TeV, collected with the CDF II detector at the Fermilab 
Tevatron. Exploiting the largest available charmed baryon sample, we measure masses and decay 
widths with uncertainties comparable to the world averages for Ec states, and significantly smaller 
uncertainties than the world averages for excited Aj states. 

PACS numbers: 14.20.Lq, 14.20.Gk 
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Hadrons containing a 5 or c quark are referred to as 
heavy-quark hadrons and provide an interesting labora- 
tory for studying and testing quantum chromodvnam- 
ics (QCD), the theory of strong interactions [l|, Be- 
cause the strong coupling constant as is large for in- 
teractions involving small momentum transfer, masses 
and decay widths of the heavy-quark states cannot be 
calculated within the framework of perturbative QCD. 
As a result, many different approaches have been devel- 
oped, for example, based on heavy-quark effective the- 
ory (HQET) 3 , nonrelativistic and relativistic potential 
models [1], or lattice QCD Q. 
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In the limit of HQET, heavy-quark mesons, comprised 
of one heavy and one hght quark, are the closest analogy 
to the hydrogen atom, which provided important tests of 
quantum electrodynamics. Heavy-quark baryons, com- 
prised of one heavy and two light quarks, extend the hy- 
drogen atom analogy of HQET by treating the two light 
quarks as a diquark system. This leads to degenerate 
spin-1/2 states resulting from the combination of a spin- 
0, or a spin-1, light diquark with the heavy quark, and 
thus represents a complementary situation compared to 
heavy-quark mesons. Measurements of the mass spec- 
trum and spin splittings of heavy-quark baryons are im- 
portant for validating the theoretical techniques, and 
build confidence in their predictions for other heavy fla- 
vor studies. 

In this paper, we measure the properties of heavy- 
quark baryons that contain a c quark, namely the 
resonances Ac(2595)+, Ac(2625)+, Ec(2455)++'°, and 
I]c(2520)'*'"'''°. For simplification, we refer to as 
Sc wherever this information is not crucial. Through- 
out the paper, the use of a specific particle state implies 
the use of the charge-conjugate state as well. The quark 
model predicts the Ac(2595)+ and Ac(2625)+, referred 
to as A*+, to be the lowest orbital excitations of the A+ 
groundstate with a spin-0 light diquark. The two Sc res- 
onances are expected to have no orbital excitation and a 
spin-1 light diquark. 

TABLE I: Theoretical predictions for the masses of the 
charmed baryons under study. All values are given in MeV/c^ . 



Hadron 


16] 


\7] 


f8, 9] 


flO] 




Ec(2455) 


2452 


2455 


2439 


2400 ± 310 


2393 


Ec(2520) 


2538 


2519 


2518 


2560 ± 240 


2489 


Ac(2595) + 




2625 


2598 


2530 ± 220 




Ac(2625)+ 




2636 


2628 


2580 ± 240 





Some theoretical predictions of the resonance masses 
are summarized in Table U where Ref. Q uses lattice 
QCD, Refs. 0-11 are based on the quark model, Ref. [l3| 
employs QCD sum rules and Ref. [ll| uses a bag model. 
There are a few calculations that predict the ILf2455) 
natural width in the region of 1-3 MeV/c^ fl^ - [l7j and 
the Sc(2520) width to be about ISMeV/c^ [ij No pre- 
dictions are available for the Ac(2595)"'" and Ac(2625)+ 
widths. 

Experimental observation of all four states studied 
here and measurements of some of their properties have 
been reported earlier [l9l - [26| . We list the world aver- 
age masses and widths in Table [Hi omitting states, 
which are difficult to reconstruct with the CDF II de- 
tector due to the inefficiency in 7r° identification. For 
I]c(2455), many measurements exist with most of the 
information coming from CLEO [H and FOCUS [13 . 
Experimental information on the Ec(2520) states comes 
exclusively from CLEO [U (2^ and it is worth noting 
that the two measurements of the Ec(2520)+"'" mass are 
inconsistent. For Ac(2595)+ and Ac(2625)+ three exper- 



TABLE II: World average values of the mass differences be- 
tween the charmed baryon resonances and the A^t mass, AM, 
and their natural widths, F [Tsf. 



Hadron 


AM [MeV/c^ ] 


F [MeV/c^ ] 


Ec(2455)++ 


167.56 ±0.11 


2.23 ±0.30 


Ec(2455)° 


167.30 ±0.11 


2.2 ±0.4 


Ec(2520)++ 


231.9 ±0.6 


14.9 ± 1.9 


Ec(2520)'' 


231.6 ±0.5 


16.1 ±2.1 


Ac (2595)+ 


308.9 ±0.6 




Ac (2625)+ 


341.7 ±0.6 


< 1.9 at 90% C.L. 



iments have contributed, namely ARGUS [2^, E687 at 
Fermilab [13, [H] and CLEO [H, ah of which suffer from 
rather small data samples. In addition, Blcchman and 
co-authors [27'| showed that a more sophisticated treat- 
ment of the mass line shape, which takes into account 
the proximity of the Ac(2595)+ mass to the sum of the 
masses of its decay products, yields a Ac(2595)+ mass 
which is 2-3MeV/c^ lower than the one observed. The 
Sc(2455) and Sc(2520) decay directly to A+tt, whereas 
the Ac(2595)+ and Ac(2625)+ end mainly in a A+tttt final 
state with dominating decays through intermediate Sc 
resonances. Therefore, these four resonances contribute 
to each other's background, which requires a dedicated 
cross-feed background modeling in each case. 

In this analysis, we exploit a large sample of A+ — > 
p K~ 7r+ decays produced in pp collisions at = 
1.96 TeV and collected by the CDF II detector. Measure- 
ments of the masses and widths of the charmed baryons 
are performed through fits to the reconstructed mass dis- 
tributions calculated from the momenta of the final state 
tracks. We take into account all expected cross-feeds and 
threshold effects. 

The paper is organized as follows. In Sec.|lT]we briefly 
describe the CDF II detector and the trigger components 
important for this analysis. Sees. IIIII and IIVI describe the 
candidate reconstruction and selection, respectively. In 
Sec. |V]we explain the fits involved in the measurements, 
followed by a discussion of systematic uncertainties in 
Sec. lVIl Finally the results and conclusions are presented 
in Sec. |VlIl 



II. CDF II DETECTOR AND TRIGGER 

Amon g th e components and capabilities of the CDF II 
detector [28] , the tracking system is the one most relevant 
to this analysis. It lies within a uniform, axial magnetic 
field of 1.4 T strength. The inner tracking volume up to a 
radius of 28 cm is comprised of 6-7 layers of double-sided 
silicon microstrip detectors An additional layer of 

single-sided silicon is mounted directly on the beam-pipe 
at a radius of 1.5 cm, allowing excellent resolution on 
the impact parameter dp i defined as the distance of clos- 
est approach of the track to the interaction point in the 
plane transverse to the beam line. The silicon detector 
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provides a vertex resolution of approximately 15 jim in 
the transverse and 70 jim in the longitudinal direction. 
The remainder of the tracking volume from a radius of 
40 to 137 cm is occupied by an open-cell drift chamber 
(COT) 30] , providing a transverse momentum resolution 
of (t(pt)/Pt ~ 0.1%/(GeV/c). Hadron identification, 
which is crucial for distinguishing slow kaons and pro- 
tons from pions, is achieved by a likelihood combination 
of information from a time-of-flight system (TOF) [sij 
and ionization energy loss in the COT. This offers about 
1.5(7 separation between kaons, or protons, and pions. 

A three-level trigger system is used for the online event 
selection. The most important device for this analysis at 
level 1 is the extremely fast tracker (XFT) [s^. It iden- 
tifies charged particles using information from the COT 
and measures their transverse momenta and azimuthal 
angles around the beam direction. The basic require- 
ment at level 1 is two charged particles with transverse 
momentum, px, greater than 2GeV/c. At level 2, the 
silicon vertex trigger [s^ adds silicon hit information 
to the XFT tracks, thus allowing the precise measure- 
ment of impact parameters of tracks. The two level 1 
tracks are required to have impact parameters between 
0.1 and 1mm and to be consistent with coming from a 
common vertex displaced from the interaction point by 
at least 100 jim in the plane transverse to the beam line. 
The level 3 trigger is implemented in software and pro- 
vides the final online selection by confirming the first two 
trigger-level decisions using a more precise reconstruction 
similar to the offline software. This trigger is designed to 
collect hadronic decays of long-lived particles such as b 
and c hadrons. As determined by a study of the impact 
parameter distributions, the sample of charmed baryons 
recorded by the trigger consists of approximately equal 
contributions from Af, decays and direct cc production. 



III. DATA SET AND RECONSTRUCTION 

The analysis is performed on a data set collected by 
the CDF II detector at the Tevatron pp collider between 
February 2002 and June 2009 corresponding to an inte- 
grated luminosity of 5.2 fb~^. The data were accumu- 
lated using the displaced two track vertex trigger de- 
scribed in the previous Section. 

The offline reconstruction of candidates starts with re- 
fitting tracks using pion, kaon and proton mass hypothe- 
ses to properly take into account differences in the mul- 
tiple scattering and ionization energy loss. In the second 
step, three tracks, one with pion, one with kaon, and one 
with proton mass hypotheses, are combined to form a A+ 
candidate. The three tracks are subjected to a kinematic 
fit that constrains them to originate from a common ver- 
tex. We require that the proton and pion candidates 
have the same charge and that the total charge of all 
three tracks is ±1. To construct Ec(2455) and Ec(2520) 
candidates we combine each A+ candidate with one of the 
remaining tracks in the event using a pion mass hypoth- 
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FIG. 1: (color online) The mass distribution of candi- 
dates used to train one of the two neural networks for the 
selection. 



TABLE III: Inputs to the neural network for the Aj selection 
sorted by their importance. 



Index Variable Index 



Variable 



LLp(p) 

(7L^^(A+) 

LLk{K) 
cos(<(A+,p)) 

do/o"do(7i") 



8 
9 
10 
11 
12 
13 



Pt{p) 
cos(<(A^,A)) 
PT(vr) 

Pt{K) 
do/cTdoip) 



esis. The Ac(2595)''" and Ac(2625)+ candidates are ob- 
tained by combining each A^ candidate with all possible 
oppositely charged track pairs taken from the remain- 
ing tracks in the event using the pion mass hypothesis 
for each of them. The tracks forming each baryon candi- 
date are subjected to a kinematic fit that constrains them 
to originate from a common point. In each step of the 
reconstruction, standard quality requirements on tracks 
and vertices are used to ensure well-measured masses and 
decay-positions. 

We use simulated events to estimate the detector mass 
resolutions of the charmed baryons studied here. The de- 
cays are simulated by means of the evtgen package [13] , 
where the A+ is forced to decay into pK ~Tr~^ with its 
resonance structure taken into account. Afterwards, the 
generated events are passed through the detector simu- 
lation and then reconstructed by the same software used 
for data. 



IV. CANDIDATE SELECTION 

The selection of the candidates is done in two steps. 
In each one we first impose some quality requirements 
to suppress the most obvious background. For the sur- 
viving candidates we use a neural network to distinguish 
signal from background. Since all final states feature a 
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FIG. 2: (color online) The mass distributions of candi- 
dates before (blue full squares) and after (red open triangles) 
requiring their neural network output to correspond to an a 
posteriori signal probability greater than 2.5%. The vertical 
dashed lines indicate a ±10MeV/c^ region around the nom- 
inal mass used for the selection of the Ec and A*"*" 
states. 



A+ daughter, the first step is the A+ selection. In the 
second step, we perform a dedicated selection of the four 
states under study. All neural networks are constructed 
with the NeuroBayes package [35, 3(| and trained, only 
using data, by means of the s'Plot technique [s^,!!!]. This 
technique assigns a weight to each candidate proportional 
to the probability that the candidate is signal. The can- 
didate weight is based on the discriminating variables, 
which are required to be independent of the ones used 
in the neural network training. In our case, the discrim- 
inating variable is the invariant mass of the candidate. 
In the training, each candidate enters with a weight cal- 
culated from the signal probability that is derived from 
its mass. Based on these weights, the neural network 
can learn the features of signal and background events. 
Since we use only data for the neural network trainings, 
we randomly split each sample into two parts (even and 
odd event numbers) and train two networks. Each of 
them is then applied to the complementary subsample 
in order to maintain a selection which is trained on a 
sample independent from the one to which it is applied. 
This approach avoids a bias of the selection originating 
from statistical fluctuations possibly learnt by the net- 
work. Additionally, using candidates from two different 
mass regions populated by background only for the train- 
ing, we verify that the network selection does not depend 
on the mass or create an artificial excess in the spectrum. 



A. At selection 



The A+ — 7> pK~TT'^ candidates are required to decay 
to a proton with pT > 1.9GeV/c and other particles 
with pt > 400MeV/c. The displacement of the associ- 
ated secondary vertex, projected onto the A^ transverse 



momentum direction, to the beam, L^y, is required to be 
greater than 0.25 mm. In addition, we use particle iden- 
tification information from the TOF and dE/dx from 
the COT. We combine the two sources of information for 
each track t into a single variable 



P. 



dE/d: 



Xi)PTOF{t) 



'l2j=Tv,K,p /j^dB/dx (^)^TOf(^) 



(1) 



where the index i denotes the hypothesis of the par- 
ticle type. The Ptof{^) is the probability to observe 
the measured time-of- flight given a particle of type i, 
and P^E/dx^^) '^^ probability to observe the mea- 
sured dE/dx. The fractions fj are /^r = 0.7, fx = 0.2, 
and Jp = 0.1, as estimated from TOF information of a 
generic background sample. We apply the requirement 
LLp > 0.6 on the proton track and LLk > 0.2 on the 
kaon track. In case TOF or dE/dx information is not 
available for a given track, we do not impose the corre- 
sponding requirement. The mass distribution of the can- 
didates with even event numbers is shown in Fig.[TJ A fit 
with a Gaussian signal and a linear background function 
defines the probability density functions (PDFs) used to 
calculate the sVlot weights for the A+ network training. 
The corresponding distribution of odd-numbered events 
is similar. 

The full list of input quantities of the neural network, 
sorted by their importance, can be found in Table Hill In 
the table, denotes the impact parameter with respect 
to the primary vertex of the pp interaction for a track 
in the plane transverse to the beam direction, its 
uncertainty, x^(Aj) the quality of the kinematic fit of the 
A+ candidate, and cos(<(A+,t)) the cosine of the angle 
between the momentum of the A^ candidate in the lab 
frame and the momentum of the proton or kaon track in 
the A+ rest frame. These angles carry information about 
the resonant substructure of the decay A+ — > pK^ 7r+. 

To demonstrate the ability of the neural network to 
classify signal and background, the mass distributions of 
A+ candidates with even event numbers before and after 
requiring their neural network output to correspond to 
an a posteriori signal probability greater than 2.5% is 
shown in Fig. [51 This requirement leads to a background 
reduction of 32% while keeping 97% of the signal. We 
use the output of the A+ neural network as input to the 
neural networks for selecting the Ec and A*"*" resonances. 



B. Ec(2455) and Sc(2520) selection 

The Ec(2455)++'" and 2^(2520)++-° ^ A+ 7r+'- se- 
lection starts with the application of a few loose require- 
ments to remove the most obvious background, followed 
by the use of a neural network. We require the a poste- 
riori signal probability of the AJ neural network to be 
greater than 2.5% (see Fig. the prij^) of the added 
pion to be greater than 400MeV/c, do{'K) < 1.5mm, and 
the mass of the A+ candidate to be within ±10MeV/c^ 
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FIG. 3: (color online) The mass difference distributions of the A^tt"*" (left) 
squares) and after (red open triangles) applying the neural network selection. 
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FIG. 4: (color online) The mass difference distribution of the 
A^TT^TT^ candidates before (blue full squares) and after (red 
open triangles) applying the neural network selection. 



well. The gVlot weights are determined by a fit to the 
AM distribution with a Gaussian function for the sig- 
nal and a linear function for the background PDF. We 
choose the threshold on the output of the Ec neural net- 
work to maximize S/\/S + B, where S is the number 
of signal Sc events and B is the number of background 
events in AM between 162.3 and 172.3 MeV/c^. The S 
and B yields are derived from a fit to the AM distribu- 
tion which uses a Gaussian function for the signal and 
a linear function for the background and covers the AM 
range used for the neural network training. The resulting 
neural network output requirement is the same for both 
charge combinations and corresponds to an a posteriori 
signal probability of the neural networks greater than 
10%. The AAf distributions of the selected candidates 
are shown in Fig. |3l 



of the nominal A+ mass Ts^, 2276.46 < M{pK- Tr+) < 
2296.46 MeV/c^ (see Fig. These requirements are 
common for both neutral and doubly-charged states. The 
mass difference AM = M{T,c) - M{A+) distributions of 
all the A+TT"*" and A+tt" candidates are shown in Fig. O 
In the AM definition, M(I]c) and M(A+) correspond to 
the reconstructed masses of the Ec and A+ candidates. 

The neural network for the final selection of the 
Sc(2455) and Ec(2520) candidates uses five input quan- 
tities. Ordered by their importance, these are the output 
of the A+ neural network NN{A'^), the proper decay 
time of the Ec candidate t{T,c) = {L^y{Tsc) ■ M(Ec))/(c • 
Pt(Ec)), the quality of the kinematic fit of the Ec candi- 
date x^(Ec), the uncertainty of the Ec impact parameter 
in the transverse plane ada{Yjc)i and the impact param- 
eter in the transverse plane of the pion from the Ec de- 
cay cJo(7r). Independent neural networks are employed 
for E++ and Ej?. The training itself is performed us- 
ing candidates in the mass difference region from 155 to 
180MeV/c2. Ahhough this includes only Ec(2455) can- 
didates, it is applied to select Ec(2520) candidates as 



C. Ac(2595)+ and Ac(2625)+ selection 

The initial step of the Ac(2595)+ and Ac(2625)+ 
A+ 7r+ TT^ selection requires the a posteriori signal prob- 
ability of the A+ neural network to be greater than 2.5%, 
2276.46 < M{pK-TT+) < 2296.46 MeV/c^ (see Fig. 
Pt{t^) of both added pions to be greater than 400MeV/c, 
and the impact parameter of the object constructed from 
the two additional pions to be (io(7r"'"7r~) < 1.0 mm. The 
mass difference AAf = A/(Ac+) - M(A+) distribution is 
shown in Fig. U) 

We use the AM region between 327 and 357MeV/c2 
for the neural network training. Although this in- 
cludes only Ac(2625)+ candidates, it is applied to se- 
lect Ac(2595)+ candidates as well. The sVlot weights 
are based on a fit to the AM distribution with a Gaus- 
sian function for the signal and a linear function for the 
background PDF. The neural network uses four inputs. 
Ordered by their importance, these are the quality of the 
A*+ kinematic fit x^lA*"*"), the uncertainty of the impact 
parameter of the combined two-pion object ado{T^^T^~), 
the output of the A+ neural network A^A''(A+), and the 
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proper decay time of the A*+ candidate t(A*+). We 
choose the requirement that maximizes S/^/S + B, cor- 
responding to an a posteriori signal probabihty of the 
neural network greater than 12.5%. The S and B yields 
are derived from a fit to the AM distribution using a 
Gaussian function for the signal and a linear function for 
the background, where we consider events in the region 
336.7 < AM < 346.7 MeV/c2. The resuhing mass differ- 
ence distribution after the final requirements is shown in 
Fig.m 



V. FIT DESCRIPTION 

To determine the mass differences relative to the 
and the widths of the six studied states, we perform 
binned maximum likelihood fits to three separate mass 
difference distributions. The first two are A^tt^ and 
A+TT", where the states £^(2455)++'° and £^(2520)++'° 
are studied. The last one is A+tt+tt^ for Ac(2595)+ and 
Ac(2625)''". In the case of the Ec states, part of the back- 
ground comes from A*"*" decays and thus has different 
properties compared to the combinatorial background. 
On the other hand, when fitting A*"*" states, there is a 
background contribution from random E++''^7r~'^ com- 
binations which have a threshold close to the Ac(2595)"'' 
state. 

The negative logarithm of the likelihood function has 
a general form of 



\nCia) = - ^In 



,7 



(2) 



where a are the free parameters, J is the number of bins 
in the histogram of the corresponding mass difference 
distribution, Uj is the number of entries in bin j, and fXj 
is the expected number of entries in bin j. The values /Xj 
are obtained using the function 

^(AM) = iVi • si(AM) + N2 ■ S2(AM) + 6(AM), (3) 

where si(AM) and S2(AM) are the PDFs for the two 
signals, fe(AM) is the background function and Ni are 
the corresponding numbers of events. All three PDFs 
depend on a subset of the free parameters a. The function 
is evaluated at the bin center to calculate the expectation 
for /ij . While the general structure is the same in all three 
fits, the PDFs are specific to Sc and A*+ states. 



A. Ec(2455) and Sc(2520) fit 

In each of the two distributions we need to parametrize 
two signals and several background components. We use 
a 150-320 MeV/c^ range to avoid complications arising 
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FIG. 5: Ec(2520)'' mass resolution obtained from sim- 
ulated events by subtracting the generated mass dif- 
ference AMgen(Ec(2520)'') from the reconstructed one 
AMroc(Sc(2520)°). The fitted function is a combination of 
three Gaussians with mean zero. 



from the description of the steep rise of the background 
at threshold. Both Ec(2455) and I]c(2520) are described 
by a nonrelativistic Breit-Wigner function. 



dN 



dAM (AM - AMo 



r2/4' 



(4) 



convolved with a resolution function. The resolution 
function itself is parametrized by three Gaussians with 
mean zero and the other parameters derived from simu- 
lated events. The average width of the resolution func- 
tion is about 1.6MeV/c2 for Ec(2455)++^" and about 
2.6MeV/c2 for Sc(2520)++^°. For illustration, the simu- 
lated I]c(2520)'^ mass resolution is shown in Fig. [S] 

We introduce a single common scaling factor s for the 
widths of all three Gaussians to correct for a possible 
mismatch in our mass resolution estimate. This scaling 
factor is allowed to fioat within a Gaussian constraint in 
the fit, what corresponds to adding 



0.5 • 



s — /i 



(5) 



with /i = 1 and a = 0.2, refiecting a 20% uncertainty on 
the mass resolution (see Sec. IVI[) . to the negative loga- 
rithm of the likelihood. 

Three different types of background are considered, 
namely, random combinations without real A+, combi- 
nations of real A+ with a random pion, and events due 
to the decay of A*"*" to A+tt+tt". The random combina- 
tions without a real Ajj" dominate and are described by 
a second-order polynomial with shape and normalization 
derived in a fit to the AM distribution from the A+ mass 
sidebands 2261.46 < M{pK-n+) < 2266.46 MeV/c^ and 
2306.46 < M{pK-n+) < 2311.46 MeV/c^. In the fit, 
this contribution is allowed to fioat within a Gaussian 
constraint implemented by the addition of 



0.5 -A^ -V- 



(6) 
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FIG. 6: Fit to the M{pK--K+ tt^ 
from mass sidebands. 



'M{pK 7r+) (left) and M(piC 7r+ tt )-M{pK 7r+) (right) distributions of the candidates 



to the negative logarithm of the likelihood, where V is 
the covariance matrix of the fit to the AAf distribution 
from the A+ mass sidebands and A is the vector of pa- 
rameters of the second-order polynomial. The fits to the 
distributions from the A^ mass sidebands can be found 
in Fig. ini The difference between doubly-charged and 
neutral spectra is due to £)*(2010)+ — > Z?" 7r+ mesons 
with multibody decays, where not all decay prod- 
ucts are reconstructed. In order to describe this reflec- 
tion, an additional Gaussian function is used. The second 
background source consisting of real A^ combined with 
a random pion is modeled by a third-order polynomial, 
where all parameters are left free in the fit. The back- 
ground originating from A*+ decays is described using 
theoretical considerations. With good approximation, 
there are two states that contribute, namely Ac(2595)+ 
and Ac(2625)"'", decaying into a A+tt+tt" final state. The 
Ac(2595)"'" decays dominantly to a ScTt final state [l^ 
and thus contributes mainly to the signal. We therefore 
neglect its contributions to the backgrounds in the AcTt 
distributions. On the other hand, the Ac(2625)+ decay is 
dominantly nonresonant 18]. To model it, we start from 
a flat A+TT+TT^ Dalitz plot and project it on the appro- 
priate axis. Since the shape of the projection depends on 
the reconstructed Ac(2625)"'" — >■ A+ 7r+ 7r~ mass, we use 
ten different values of this mass and weight their contri- 
bution according to the Ac(2625)+ shape we obtain from 
our flt to the A+tt+tt" data. This contribution amounts 



to about 2% of the total background. 

The full flt to the AAf distribution, containing all sig- 
nal and background components, can be found in Fig. [T] 
The value of the E++ flt is 340 (324 degrees of free- 
dom) and that of the I]° fit is 384 (321 degrees of free- 
dom). 



B. Ac(2595)+ and Ac (2625)+ flt 

The fit for Ac(2595)+ and Ac(2625)+ includes two 
signals and several background components and is per- 
formed in a AM region from 290 to 400MeV/c^. Pre- 
vious measurements of the Ac(2595)+ properties indi- 
cate that it decays dominantly to the final state EgTr, 
with the threshold very close to the Ac(2595)+ mass p^ . 
This introduces an additional complication to the fit com- 
pared to the Ec case. Blechman et al. [275 showed that 
taking into account the mass dependence of the natu- 
ral width yields a lower Ac(2595)+ mass measurement 
than observed by previous experiments. With the present 
event sample we are more sensitive to the details of the 
Ac(2595)+ line shape than previous analyses and include 
this dependence in the model. 

The Ac(2595)+ parametrization follows Ref. '27]. The 
state is described by a nonrelativistic Breit-Wigner func- 
tion of the form 



J 



dN r(A+7r+7r-) 
(X - - - 

dAM (AAf - AMa^(2595)+)' + (r(A+^+7r-) + r(A+7r"^0))2/4 ' 

I 



(7) 



where r(A+7r+7r^) and r(A+7r"Tr") are the mass- 
dependent partial widths to the A+tt+tt" and A+7r°7r° 
final states. Assuming that those two final states satu- 



rate nearly 100% of the Ac(2595)"'" decay width, the sum 
in the denominator corresponds to the total width. The 
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FIG. 7: (color online) The M{pK~ tv+ tv+) - M{pK~-K+) (left) and M{pK~n+ Tr~) - M{pK~-K+) (right) distributions obtained 
from data (points with error bars) together with the fits (black solid line). The brown dashed and purple dotted lines correspond 
to the two signal contributions, the green dash-double-dotted line represents the combinatorial background without real A+, 
the blue long-dashed line shows real combined with a random pion and the red long-dash-dotted line represents a reflection 
from A*+ decays. The red dash-dotted line corresponds to the sum of all three background contributions. 



two partial widths are derived in Ref. (39| as 

r(A+7r+7r-) = ^g||^m^+ J dE,dE2{\p2f\A{E,)f 
+ |pi|2|B(S2)|' + 2pi • p2Re[AiEi)B*iE2)]), 



(8) 



5l 



16^3/4 



'At 



dE,dE2i\p2\^\CiEi 



+ |pinc(^2)|' + 2pi • p2MCiEi)C*iE2)]). 

(9) 

Here, = 132MeV/c^ is the pion decay constant (40| . 

is the world average A+ mass, Ei,E2 are the ener- 
gies of the two pions in the rest frame of the Ac(2595)^, 
and Pi,p 2 are the corresponding momenta. Following 
Ref. [27| . the coupling constant 52 is determined by the 
Sc decay width using the relation 



5i 



27r/2 ms. 



\Pn\ 



(10) 



with TOs^ being the world average mass of the I]c(2455) 
and the momentum of the pion from the I]c(2455) 
decay to AcTt in the Ec(2455) rest frame. From the 
world average Fs^ — 2.2MeV/c^ [l8| we obtain the value 
§2 = 0.365 which is fixed in the fit. The amplitudes A, 
B, and C for the decays Ac(2595)+ ^ 2^(2455)° 7r+, 
Ac(2595)+ ^ Sc(2455)++7r-, and Ac(2595)+ ^ 
Sc(2455)+7r" are parametrized as 

h2E 



A{E) 



Am - Am^o - E + iF^o /2 ' 



(11) 



B{E) 



h2E 



Am — Atoj2++ — E + /2 



1 



CiE) = - 



hoE 



2 Am - Amj,+ - E + iT^+ /2 ' 



(12) 



(13) 



In these definitions, m-^++,+,o and Fj^++,+,o are the mass 
and the width of the Ec(2455)++'+'° taken from Ref. [H]. 
The coupling constant /12, defined in Ref. [l^, is related 
to the decay width of the Ac(2595)''" and represents the 
actual quantity we measure instead of the natural width. 
This approach describes a purely S'-wave decay, a pos- 
sible ZJ-wave contribution is assumed to be negligible 
and ignored. For illustration, we show the dependence 
of the two partial widths on M(Ac(2595)+) - M(A+) 
in Fig. [SI The shape defined by Eq. [7] is then numer- 
ically convolved with a resolution function determined 
from simulation and consisting of three Gaussians with 
mean zero. The average width of the resolution function 
is about 1.8MeV/c^. As for the Ec case, we introduce 
a common, Gaussian constrained, scaling factor for the 
widths of all three Gaussians, in order to account for the 
uncertainty in the width of the resolution function. 

The signal PDF for the Ac(2625)+ is the nonrelativis- 
tic Breit-Wigner function of Eq.|4] convolved with a three 
Gaussian resolution function determined from simula- 
tion, which has an average width of about 2.4MeV/c^. 
Again, all three Gaussians have mean zero and a com- 
mon, Gaussian constrained, scaling factor for their 
widths is introduced. 

The background consists of three different sources, 
which include combinatorial background without real 
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FIG. 8: (color online) Calculated dependence of r{A'^TT'^Tv ) 
and r(A+7rV°) on M(Ac(2595)+) - M(A+). The constant 
factor ^2 is determined by a fit to the experimental data. 




FIG. 9: Fit to the M(piC"7r+ tt+tv') - M(pA'"7r+) distribu- 
tion of the candidates from Ajf mass sidebands. 

A+, real combined with two random pions, and real 
S^^'" combined with a random pion. The combina- 
torial background without real Ajj" is parametrized by 
a second-order polynomial whose parameters are deter- 
mined in a fit to the AM distribution of candidates 
from the A+ mass sidebands, 2261.46 < M{pK~TT+) < 
2266.46 MeV/c2 and 2306.46 < M{pK-TT+) < 
2311.46 MeV/c^. This distribution is shown in Fig. M 
together with the fit. In the final fit, we keep the param- 
eters for this background fioating within a Gaussian con- 
straint of the form of Eq.[6]to the values found in the fit to 
the candidates from the A+ mass sidebands. The second 
source, consisting of real A^ combined with two random 
pions, is parametrized by a second-order polynomial with 
all parameters allowed to float in the fit. The final source 
of background are real Sc combined with a random pion. 
For this source, the main issue is to have the proper shape 
close to the threshold. Small imperfections at higher 
AM can be ignored, as the second background source 
has enough flexibility to absorb it. The PDF of this Ec 
background is based on a constant function defined from 
the threshold to the end of the fit range. In order to take 



into account the natural widths as well as resolution ef- 
fects, we use the weighted sum of ten such functions for 
both 5]c(2455)++ and Ec(2455)°. Their thresholds and 
weights are chosen according to the shapes derived in the 
Ec fits shown in Fig. [T] The size of this contribution is 
constrained to the Ec(2455) yield obtained from the fits 
to the M(Ec) — M{pK~7T~^) distributions for candidates 
with M(A*+) - M{pK-Tr+) > 355MeV/c2. These two 
distributions together with the fits are shown in Fig. [TO] 
The full fit to the AM distribution, containing all sig- 
nal and background components, can be found in Fig. II II 
The value of the fit is 227 (206 degrees of freedom). 
Compared to that, the value of a fit with a mass- 
independent Ac (2595)+ decay width, shown in Fig. [121 
increases to 286 (206 degrees of freedom). 

VI. SYSTEMATIC UNCERTAINTIES 

We investigate several systematic effects that can affect 
the measurements. Generally, they can be categorized as 
imperfect modeling by the simulation, imperfect knowl- 
edge of the momentum scale of the detector, ambiguities 
in the fit model, and uncertainties on the external inputs 
to the fit. In this Section we discuss how they can affect 
our results and the way we assess them. A summary of 
the assigned uncertainties can be found in Tables HVlWII 
To obtain the total systematic uncertainties, we add up 
the contributions from all sources in quadrature. 



A. Mass resolution model 

To properly describe the signal shapes, we need to 
understand the intrinsic mass resolution of the detec- 
tor. Since we estimate this using simulated events, 
it is necessary to verify that the resolution obtained 
from simulation agrees with that in real data. We use 
£)*(2010)+ D°t:+ with D° K' tt+ decays and 
'0(25') J/i/'TT+TT^ with J/0 /-^^ decays for this 

purpose. We compare the resolution in data and simu- 
lated events as a function of the px of the pions added 
to 13° or J/ip as well as the instantaneous luminosity. We 
also compare the overall resolution scale between data 
and simulated events and find that all discrepancies are 
less than 20%, which we assign as uncertainty on our 
knowledge of the resolution function. The contribution 
from this uncertainty is already included in the uncer- 
tainties on the resonance parameters determined by the 
default fit with Gaussian constraint on the resolution 
scaling factor s, the resulting values for which are listed in 
Table |Vni These values are consistent with 1, indicating 
that the resolution is well understood within the assigned 
uncertainty. To disentangle it from the statistical com- 
ponent, we repeat the fits on data without multiplying 
the widths of the resolution function by the scaling fac- 
tor s from Eq. [5] The systematic uncertainty due to the 
imperfect modeling of the resolution function is then ob- 
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TABLE IV: Systematic uncertainties on the measurements of the mass diflerences and decay widths of the E^"'' resonances. 
The corresponding statistical uncertainties are listed for comparison. 



Source 


AM(Ec(2455)++) 
[MeV/c^ ] 


r(Ee(2455)++) 
[MeV/c^ ] 


AM(Ec(2520)++) 
[MeV/c^ ] 


r(Ee(2520)++) 
[MeV/c^ ] 


Resolution model 




0.40 




0.69 


Momentum scale 


0.12 


0.20 


0.12 


0.20 


Fit model 


0.02 




0.11 


1.16 


External inputs 










Sum 


0.12 


0.45 


0.16 


1.36 


Statistical 


0.04 


0.13 


0.56 


2.12 



tained by the difference in quadrature of the uncertainty 
of the fit with and without the Gaussian constraint. This 
uncertainty in the resolution has a large impact on the 
natural widths, but a negligible effect on the mass differ- 
ences. 



B. Momentum scale 

The accuracy of the momentum scale depends on the 
precision with which the magnetic field and the amount 
of material in the detector are known. Both effects are 
originally calibrated using J/i/; A*^ decays [4l[. A 

limitation of this calibration is that it uses muons that 
are required by the detector acceptance to have pT > 
1.5GeV/c, while pions from Sc or A*+ decays typically 
have much lower px- The estimate of the uncertainty on 
the mass differences comes from our previous work on the 
X(3872) hadron There, i}{2S) J/iPn+n- decays 
are used to study the momentum scale uncertainties by 
comparing the measured 'ijj{2S) mass with the world aver- 
age value pL^. In addition, we study the ip{2S) mass de- 
pendence on the kinematic properties of the pions, which 
constrains the sizes of possible effects. Furthermore, we 
verify the momentum scale by using D* (2010)+ — > tt"*" 
decays, where the resulting deviation from the world av- 
erage is far below the uncertainty derived from ■0(25'). 
Based on Ref. [42|, we assign a 0.12MeV/c^ uncertainty 



on the mass differences of all states under study due to 
the imperfect knowledge of the momentum scale. The 
corresponding effect on the natural widths was studied 
in our previous measurements of the masses and widths 
of the excited charmed meson states and ' 
we assign the 0.2MeV/c^ found there as the uncertainty 
on the natural widths due to this source. To translate 
this uncertainty to the coupling constant /12, we assign 
it to the sum r(A+7r+7r") + r(A+7r°7r°) (see Eqs. |5]and 
|9]) , which is a function of h2 , and perform Gaussian error 
propagation. 



C. Fit model 

In terms of our fit model and procedure we check two 
effects, the internal consistency of the fit and the shape of 
the signal PDFs. We do not perform an explicit check of 
the background parametrizations as those are described 
by polynomials and any analytic function can be approx- 
imated by a polynomial of sufficient complexity. Since 
the fit quality does not indicate significant discrepancies 
between data and the model, we conclude that the degree 
of the polynomial functions used is sufficient. Some back- 
grounds are determined from independent sources, but as 
the appropriate parameters are Gaussian constrained in 
the fit, the uncertainty originating from the sample size of 
the external sources, like A^ mass sidebands, is already 
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TABLE V: Systematic uncertainties on the measurements of the mass differences and decay widths of the E^! resonances. The 
corresponding statistical uncertainties are hsted for comparison. 



Source 


AM(2jc(2450j ) 

[MeV/c^ ] 


1 (Lc(2455) ) 
[MeV/c^ ] 


AM(Ec(2520)") r(Ec(2520)") 
[MeV/c^ ] [MeV/c^ ] 


Resolution model 




0.45 




0.70 


Momentum scale 


0.12 


0.20 


0.12 


0.20 


Fit model 


0.02 




0.11 


1.16 


External inputs 










Sum 


0.12 


0.49 


0.16 


1.37 


Statistical 


0.03 


0.11 


0.43 


1.82 



TABLE VL Systematic uncertainties on the measurements of the mass differences of the A*^ resonances and the pion coupling 
constant /12 (r(Ac(2595)"'")). The corresponding statistical uncertainties are listed for comparison. 



Source 


AM(Ac(2595)+ 
[MeV/c^ ] 


) hi 


r(Ac(2595)+ 
[MeV/c^ ] 


) AA/(Ac(2625)+) 
[McV/c^ ] 


Resolution model 


0.06 


0.03 


0.22 




Momentum scale 


0.12 


0.03 


0.20 


0.12 


Fit model 










External inputs 


0.15 


0.06 


0.36 




Sum 


0.20 


0.07 


0.47 


0.12 


Statistical 


0.14 


0.04 


0.30 


0.04 



TABLE VIL Mass resolution scaling factors s floating within 
Gaussian constraints in the fits. 



Hadron s 
Ec(2455)++ 0.93 ±0.17 
Ec(2455)° 1.07 ±0.13 
Ec(2520)++ 1.02 ±0.20 
Ec(2520)° 1.00 ±0.20 
Ae(2595)+ 0.95 ±0.15 



included in the statistical uncertainties of the results. 

To check the internal consistency of the fit procedure, 
we generate a large ensemble of statistical trials using 
PDFs of our fit model with parameters obtained from the 
fit to data. Estimates of all physics parameters except 
the mass differences and natural widths of the Sc(2520) 
resonances are found to be unbiased. The Ec(2520) mass 
differences have small biases towards higher values and 
the I]c(2520) natural widths are biased towards lower 
values. These biases on the I]c(2520) resonance parame- 
ters result from the fairly low signal to background ratio 
and the flexibility in the background PDF, which tends 
to absorb the tails of the relatively broad signal struc- 
ture. We repeat the study with a true value for the 
i;c(2520) natural width below (T = T.SMeV/c^) and 
above (F = 20 MeV/c^) the measured value and find that 
the biases have a small dependence on the true value. 
The biases are largest for a true value of the natural 
width of 20 MeV/c^ and we consequently assign these bi- 
ases as systematic uncertainties on the mass differences 
and natural widths of the Sc(2520) states. 

Concerning the uncertainty on the signal shape, we 
check whether our signal parametrization using nonrela- 



tivistic Breit-Wigner functions provides a proper descrip- 
tion. We refit the Ec and Ac(2625)+ data using a P-wave 
relativistic Breit-Wigner function of the form 



diV 

dm '^'^^ 



m ■ F(m) 



i2)2 .r2(m) 



with 



F(to) 



° \qo) V mJ \l + qm^ 



(14) 



(15) 



where m — AM + , R is the Blatt-Wcisskopf radius 
set to 3(GeV/c)~^ 0, Hi], mo and Fo are the nomi- 
nal mass and width, and q{qo) is the momentum of the 
daughters in the Ec or Ac(2625)"'' rest frame calculated 
from the nominal mass. For the Ac (2595)+ we replace 
the nonrelativistic Breit-Wigner function of Eq. [7] by a 
relativistic one and use the variable width defined in 
Eqs.[5]and[ni For the Ec(2455) we observe a difference of 
0.02MeV/c2 in the mass difference, which we assign as 
a systematic uncertainty. In the cases of Ec(2520) and 
A*"'' resonances we do not observe any shift and conclude 
that the effect is negligible. 



D. External inputs 

Finally, the line shape of the Ac(2595)"'" depends on the 
input values of the Ec(2455) masses and widths and the 
pion decay constant /^r. We repeat the fit using values of 
those parameters smaller or larger by 1 standard devia- 
tion and take the stronger variation as systematic uncer- 
tainty. The effect of the uncertainty on the world average 
Ec(2455) masses and widths used as input is dominant 
compared to the effect of the uncertainty on f^. 



14 



1000 



> 



o 



(0 

■g 

T3 
c 
as 
O 



800 



600 



400 



200 



- — Data 

— Fit Function 

■■■■ Aj(2595)*->A*7l*ll" 

- A ^2625)* -> a; 71* It" 




N(Ac(2595)*) = 3500 
N(A„(2625)*) = 6200 











TABLE VIII: Measured resonance parameters, where the first 
uncertainty is statistical and the second is systematic. 
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Mass(pKV ji+n")-Mass(pKV) [MeV/c ] 

FIG. 11: (color online) The M{pK- n+ n+n-) - M{pK- tt+) 
distribution obtained from data (points with error bars) to- 
gether with the fit (black solid line). The brown dashed and 
purple dotted lines correspond to the two signal contributions, 
the green dash-double-dotted line represents the combinato- 
rial background without real A^, the blue long-dashed line 
shows real A^ combined with two random pions and the red 
long-dash-dotted line represents real Ec combined with a ran- 
dom pion. The red dash-dotted line corresponds to the sum 
of all three background contributions. 
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FIG. 12: (color online) The M{pK-n+ n+Tv-) - M{pK-n+) 
distribution obtained from data (points with error bars) to- 
gether with the fit (black solid line), where a Breit-Wigner 
function with a mass-independent decay width is used to 
model the Ac(2595)^ line shape. Explanations of the vari- 
ous background contributions can be found in the caption of 
Fig. El 



Hadron 



E,(2455)++ 
Ec(2455)° 
Ec(2520)++ 
Ec(2520)° 
Ac(2595)+ 
Ae(2625)+ 



AM [McV/c^ 



F [MeV/c^ 



167.44 ±0.04 ±0.12 
167.28 ±0.03 ±0.12 
230.73 ±0.56 ±0.16 
232.88 ±0.43 ±0.16 
305.79 ±0.14 ±0.20 
341.65 ±0.04 ±0.12 



2.34 ±0.13 ±0.45 
1.65 ±0.11 ±0.49 
15.03 ± 2.12 ± 1.36 
12.51 ± 1.82 ± 1.37 
0.36 ±0.04 ±0.07 



VII. RESULTS AND CONCLUSIONS 

We perform fits to the M(pi^"7r+ 7r+) - M{pK-7r+), 
M{pK-TT+ TT-) - M{pK-Tr+), and M{pK-TT+ tt+tt") - 
M{pK^TT^) mass difference distributions to obtain the 
desired resonance properties. The data distributions and 
fits are shown in Figs. [7] and [TTJ We select about 13800 
Ec(2455)++, 15900 £^(2455)°, 8800 £^(2520)++, 9000 
Sc(2520)O, 3500 Ac(2595)+, and 6200 Ac(2625)+ signal 
events. The resonance parameters obtained can be found 
in Table rvml For the width of the Ac(2625)+ we observe 
a value consistent with zero and therefore calculate an 
upper limit using a Bayesian approach with a uniform 
prior restricted to positive values. At the 90% credibility 
level we obtain r(Ac(2625)+) < 0.97MeV/c2. For easier 
comparison to previous results [1^ H^, corresponds 
to a Ac(2595)+ decay width of r(Ac(2595)+) = 2.59 ± 
0.30±0.47MeV/c2, calculated at AM(Ac(2595)+). Our 
precise measurement of the coupling constant /12 can, for 
instance, be used to predict the width of the 2^(2645), 
as discussed in Ref. ^] . 

In Figs. I13H151 our results are compared to previ- 
ous measurements by other experiments. Except for 
AM(Ac(2595)+), all our measurements agree with the 
previous world average values. For AAf (Ac(2595)+) we 
show that a mass- independent natural width does not de- 
scribe the data (see Fig. [12]) and observe a value which is 
3.1 MeV/c^ smaller than the existing world average. This 
difference is the same size as estimated in Ref. [231 ■ Since 
this data sample is 25 times larger than the ones studied 
so far, our results on the properties of A*+ states pro- 
vide a significant improvement in precision compared to 
previous measurements. The precision for the Sc states 
is comparable to the precision of the world averages. 
Concernin g th e inconsistency of the two CLEO measure- 
ments [21!. [2^ of the I]c(2520)'^'^ mass, our data favor a 
smaller value. 

In conclusion, we exploit the world largest samples of 
excited charmed baryons to measure the resonance pa- 
rameters of six states, namely Sc(2455)+"'", Ec(2455)'', 
Sc(2520)++, Ec(2520)'', Ac(2595)+, and Ac(2625)+. Ta- 
ble IIXI summarizes the results for their masses and 
widths. These measurements provide a significant im- 
provement in the knowledge of the resonance parameters 
of the states and represent the first analysis of charmed 
baryons at a hadron collider. 
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TABLE IX: Results for the masses and widths of the charmed 
baryons under study. The first uncertainty is the combined 
statistical and systematic experimental uncertainty. For the 
masses, the second uncertainty originates from the world av- 
erage Ajt mass 18]. 



Hadron 



Ec(2455)++ 
Ec(2455)° 

Ec(2520)++ 
Ec(2520)° 
Ac(2595)+ 



M [MeV/c- 



12- 



r [MeV/c^ ] 



2453.90 ±0.13 ±0.14 
2453.74 ±0.12 ±0.14 
2517.19 ±0.46 ±0.14 
2519.34 ±0.58 ±0.14 
2592.25 ±0.24 ±0.14 



2.34 ±0.47 
1.65 ±0.50 
15.03 ±2.52 
12.51 ± 2.28 



Ac(2625)+ 2628.11 ± 0.13 ± 0.14 < 0.97 at 90% C.L. 
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FIG. 13: Comparison of our results for the Se(2455) mass differences and decay widths with previous measurements by 
Fermilab E791 [13, FOCUS [131, and CLEO [Q. The error bars represent the statistical (vertical marks) as well as the 
combined statistical and systematic uncertainties. The vertical dashed line together with the surrounding box symbolizes the 
world average value and its uncertainty taken from Ref. [1^. This average does not take into account the measurement at 
hand. 
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FIG. 14: Comparison of our results for the Ec(2520) mass differences and decay widths with previous measurements by 
CLEO [21I. [2^. Further explanations can be found in the caption of Fig. 1131 
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FIG. 15: Comparison of our results for the Ac(2595)"'" mass difference and decay width as well as the Ac(2625)"'" mass difference 
with previous measurements by CLEO [26[, Fermilab E687 [2^. [25l|. and ARGUS fSS]. Further explanations can be found in 
the caption of Fig. [131 



